Open heavy flavor and J/psi at RHIC and LHC 
within a transport model 



Jan Uphoff, 1 Oliver Fochler, 1 Zhe Xu, 2 and Carsten Greiner 1 

1 Institut fur Theoretische Physik, Johann Wolfgang Goethe-Universitat Frankfurt, 
Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany 

2 Department of Physics, Tsinghua University, Beijing 100084, China 

E-mail: uphof f @th . physik . uni-f rankf urt . de 

Abstract. The production and space-time evolution of heavy quarks and J/ip in the 
quark gluon plasma is studied within the partonic transport model Boltzmann Approach to 
MultiParton Scatterings (BAMPS). This framework allows interactions among all partons: 
gluons, light quarks and heavy quarks. Heavy quarks, in particular, interact with the rest of 
the medium via binary scatterings with a running coupling and a more precise Debye screening 
which is derived from hard thermal loop calculations. We compare our results of the elliptic 
flow and nuclear modification factor not only to experimental data of heavy flavor electrons at 
RHIC, but also to LHC data of heavy flavor electrons, muons, D mesons, and non-prompt J/ip. 
Where no data is available yet, we make predictions for those observables. Furthermore, results 
on prompt J/ip elliptic flow are reported for RHIC energy within the same framework, taking 
the dissociation as well as regeneration of J/tp in the quark-gluon plasma into account. 



1. Introduction 

In ultra-relativistic heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) pQ [2] and 
the Large Hadron Collider (LHC) |3j a unique state of matter is produced, in which quarks 
and gluons form the relevant degrees of freedom. This quark gluon plasma (QGP) has exciting 
properties such as collective behavior like a nearly perfect liquid or the quenching of high energy 
particles. 

In particular heavy quarks (charm and bottom) provide an insightful way to learn more about 
the properties of this matter. Since they are heavy, their production time is at a very early 
stage of the heavy ion collision when enough energy is available In consequent interactions 
with other particles they are influenced by the medium without fully thermalizing |5j and, 
therefore, carry information about its properties. Charm and bottom quarks travel through 
the medium, lose energy and participate in the collective motion. The exact mechanism of 
the intense interaction with the medium is actively debated, most recently, for instance, in 
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Complimentary to these open heavy flavor particles are particles in which the flavor is hidden. 
That is, bound states of a heavy quark and an anti-heavy quark. They can survive in the QGP 
to some extent, but melt if the temperature of the medium is too high [16]. Therefore, they 
could be used as a thermometer. In addition to this melting they can also be regenerated by 
two heavy quarks which meet in the medium. The most prominent hidden heavy flavor particle 
is the J/ip, which consists of a charm and an anti-charm quark. 



In this paper we employ the partonic transport model BAMPS to study the suppression and 
elliptic flow of both open and hidden heavy flavor. In detail, we make calculations for open heavy 
flavor observables such as the nuclear modification factor Raa and elliptic flow V2 of heavy flavor 
muons and electrons, D mesons and non-prompt J/ip which stem from B mesons. Furthermore 
we present V2 calculations of prompt J/tp. 

2. Partonic transport model BAMPS 

For the simulation of the QGP we use the partonic transport model Boltzmann Approach to 
MultiParton Scatterings (BAMPS) |17l I18j . which describes the full space-time evolution of the 
QGP by solving the Boltzmann equation, 
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for on-shell partons and pQCD interactions. On the heavy flavor sector all possible binary 
collisions are implemented, such as gg <-> QQ, qq «-> QQ, gQ — > gQ, qQ — > qQ, and 
J/vp + g-v^c + c. As a note, heavy quark production via gluon fusion within BAMPS is 
discussed in Ref. [3], but will not be addressed in this paper. 

The divergent t channel of elastic processes is regularized with a screening mass fi which 
is determined by matching energy loss calculations within the hard thermal loop approach 
|19 t I20|, [2T| . Furthermore, the running of the coupling is explicitly taken into account. We note 
that currently for light partons neither a running coupling nor an improved Debye screening 
is employed, which we plan, however, to do in the future. Inelastic processes 2 — >• 3 are 
currently being studied for heavy quarks, but have already been implemented for light partons. 
Details can be found in the following references: the model itself [17\ 118] . the open heavy flavor 
implementation @J EQ [8] , and the J/if> details [221I23] . 

3. Results and comparison with data 

The implemented running coupling and the improved screening procedure, which reproduces 
the energy loss from hard thermal loop calculations, effectively enhances the heavy quark cross 
section with the medium. Quantitative comparisons [Ml ED ESI [8] show that elastic processes 
contribute significantly to the energy loss of heavy quarks. However, they alone can reproduce 
neither the data of the nuclear modification factor nor the elliptic flow of any heavy flavor particle 
species. This is not too surprising since we expect that radiative 2—^3 processes also play an 
important role and that both processes together should account for the measured suppression 
and flow. The quantitative contribution of radiative processes will be studied in an forthcoming 
investigation. 

In this paper we mimic the radiative influence by effectively increasing the elastic cross section 
by a factor K = 3.5 which fits the heavy flavor electron elliptic flow data from PHENIX in the 
centrality class 20-40 % [8] . On the left hand side of Fig. [1] we show our prediction for 0-60 % 
centrality class for the same parameters, which is compared to new data from STAR [26J. Not 
surprisingly the agreement is very well for small and intermediate transverse momenta. Even the 
bump around 1.7 GeV is nicely described. At large pr, however, the data is significantly larger 
than our calculations. This could also be due to jet-like correlations which are not considered in 
the analysis [26]. Simultaneously, the Raa data at RHIC can also be described with the same 
parameter [8]. Whether this effective description is valid and radiative contributions really boil 
down to simply multiplying the elastic cross section by a constant factor will be studied with 
BAMPS in the future. 

On the right hand side of Fig. [1] the nuclear modification factor of D mesons for very central 
events at LHC is depicted. Our prediction with BAMPS is slightly smaller than the experimental 
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Figure 1. Left: Elliptic flow vi of heavy flavor electrons at RHIC for 0-60 % centrality class as 
a function of transverse momentum with STAR data [26]. Right: Nuclear modification factor 
Raa of D mesons for b = 3.6 fm with data [27j . 
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Figure 2. Left: Elliptic flow v% of D mesons, non-prompt J/ip, muons, and electrons at LHC 
together with available data for D mesons and electrons |28t 129] . Right: Elliptic flow of J/ip 
with data [30]. For comparison the charm quark V2 is also shown. 



data points. This is in agreement with the observation that we also underestimate the Raa of 
non-prompt J/ip, heavy flavor electrons, and muons [5]. 

The left hand side of Fig. [2] shows our predictions of the elliptic flow of D mesons, non- 
prompt J/ip, heavy flavor electrons, and muons. For D mesons and electrons there is already 
data available which agrees very well with our calculations. 

Experimental measurements at RHIC showed that the vi of prompt J/ip is very small [30] . 
This is in contradiction with the regeneration picture where the flow of the charm quarks should 
be transferred to the J/ip- BAMPS is an ideal framework to study this in more detail since it 
reproduces the open heavy flavor flow and also allows recombination of charm quarks to J/ip. 
The right panel of Fig. [2] shows that in BAMPS even with regeneration the elliptic flow of all 
J/ip is compatible with the data. 

4. Summary 

We presented heavy flavor calculations with the parton cascade BAMPS. So far only elastic 
interactions with a running coupling and an improved Debye screening are implemented on 



the heavy flavor sector. To account for missing radiative contributions we scale the elastic 
cross section with K = 3.5 which is determined from comparing to the PHENIX vi data at 
RHIC. With this parameter we find also a good agreement to recent Vi data from STAR in 
another centrality class. The predictions for all open heavy flavor Raa at the LHC consistently 
underestimate the data slightly. However, the predicted V2 of D mesons and electrons is in very 
good agreement with recent data from ALICE. The radiative processes [31] are currently being 
implemented in BAMPS. It will be interesting to see if they can indeed account for scaling 
the binary interactions with a constant factor. Furthermore, we showed some first BAMPS 
calculations on prompt J ftp vi at RHIC. 
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